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Abstract 

Neutrino oscillations are treated from the point of view of relativistic first quantized theories 

o i 

and compared to second quantized treatments. Within first quantized theories, general oscillation 

5h 

Oh" probabilities can be found for Dirac fermions and charged spin bosons. A clear modification in 

<" 

the oscillation formulas can be obtained and its origin is elucidated and confirmed to be inevitable 

from completeness and causality requirements. The left-handed nature of created and detected 

>' 

' neutrinos can also be implemented in the first quantized Dirac theory in presence of mixing; the 

o ' 

probability loss due to the changing of initially left-handed neutrinos to the undetected right- 
handed neutrinos can be obtained in analytic form. Concerning second quantized approaches, it 

in ■ 

is shown in a calculation using virtual neutrino propagation that both neutrinos and antineutrinos 



Oh 

CD 



O may also contribute as intermediate particles. The sign of the contributing neutrino energy may 

have to be chosen explicitly without being automatic in the formalism. At last, a simple second 
quantized description of the flavor oscillation phenomenon is devised. In this description there is no 
interference terms between positive and negative components, but it still gives simple normalized 
oscillation probabilities. A new effect appearing in this context is an inevitable but tiny violation 
of the initial flavor of neutrinos. The probability loss due to the conversion of left-handed neutrinos 
to right-handed neutrinos is also presented. 
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I. INTRODUCTION 

Compelling experimental evidences [j| have shown that neutrinos undergo flavor oscil- 
lations in vacuum. Consequently, this fact requires massive neutrinos with mixing. These 
ingredients are not present in the standard model of elementary particles. For this reason, 
on the one hand, neutrino oscillations can provide a direct window to probe physics beyond 
the standard model Q. On the other hand, some theoretical studies of mixing in the context 
of quantum field theory (QFT) by Blasone and Vitiello (BV) ^,0] show the mixing problem 
may be related to more fundamental issues such as unitarily inequivalent representations 
and the vacuum structure, and its study is theoretically interesting for its own sake. 

Nevertheless, the simpler plane-wave quantum mechanical descriptions |f| seemed to 
be in accordance, in certain realistic limits, with more refined descriptions, including vari- 
ous ingredients, such as localization aspects , flavor current densities [lOf, influence 
of creation and detection processes 11 , [l2 L time-dependent perturbation theory jl^ , and 
intermediate neutrinos with path integrals 14] . Moreover, many treatmens within the quan- 



tum field theory (QFT) framework were also proposed |l5l. Ilfil Il7l Il8l 1 13. 12(1 l2l|. aiming to 
solve the various unclear aspects of the quantum mechanics of neutrino oscillations jlfll . . 

It has been known for a long time that the coherence necessary for neutrino oscillations 
depends crucially on localization aspects of the particles involved in t he p roduction of neu- 
trinos 0. This point of view can be supported by QFT arguments jl£| as well. It raises 
then the question of how the coherent superposition of mass eigenstate neutrinos, which is 
called a "flavor" eigenstate, is created |2i|. One way that became customary to avoid the 
ambiguities involving the question on how neutrinos are created and detected is to use an ex- 
ternal (E) wave packet (WP) approach |l7|, in contrast to an intermediate (I) WP approach. 
According to Ref. the IWP treatments are the simpler first quantized ones treating the 
propagation of neutrinos as free localized wave packets. In contrast, EWP approaches con- 
sider localized wave packets for the sources and detection particles while the neutrinos were 
considered intermediate virtual particles. The central issue distinguishing the general IWP 
and EWP approaches is: despite its direct unobservability, is the intermediate neutrino a 
real (on-shell) particle propagating freely? If the answer is affirmative, the IWP approaches 
would be a good approximation of the oscillation phenomena. 

On the other hand, another classification scheme can be used do classify the various 
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existing treatments considering a more physical criterion irrespective of the use of WPs. 
Consider the descriptions of neutrino oscillations that (A) include explicitly the interactions 
responsible for the mixing and those (B) that only treat the propagation of neutrinos, i.e,. 
the mixing is an ad hoc ingredient. A more subtle aspect in between would be the (explicit or 
phenomenologically modelled) consideration of the production (and detection) process(es). 
In general, the IWP approaches are of type (B). The EWP approaches are of type (A). The 
BV approach, although in the QFT formalism, is of type (B) since mixing is introduced 
without explicitly including the interaction responsible for it. The type (B) approaches 
have the virtue that they can be formulated in a way in which total oscillation probability 
in time is always conserved and normalized to one 0, |^. This feature will be present 
in all first quantized approaches treated here (sees. HU and IIII|) and in a second quantized 
version f sec. IIV AT) . If different observables are considered, or a modeling of the details of the 
production and detection processes is attended, further no—on is necessary QflQ- 
In such cases, the oscillating observable might differ from the oscillation probability. On the 
other hand, type (A) approaches tend to be more realistic and can account for the production 
and detection processes giving experimentally observable oscillation probabilities j^jj. Of 
course, they are essential to the investigation of how neutrinos are produced and detected [ill 
I20I ]. We are not directly interested in these matters here. 

Considering first quantized type (B) app roaches, some recent works treating the flavor 
oscillation for spin one-half particles 22, 2^ have already find additional oscillatorial effects 
compared to usual oscillation formulas with WPs 0,0]. These effects are investigated and 
it is shown in sec. [H] how these additional oscillatorial behavior, which have characteristic 
frequencies much greater than usual oscillation frequencies, comes from the interference be- 
tween positive and negative frequency components of the initial WP. It can be understood 
as a consequence of the impossibility to simultaneously exclude all negative energy contri- 
butions of the initial spinorial wave function with respect to bases characterized by different 
masses. Moreover, this rapid oscillations are always present, independently of the initial 
WP, if a well defined flavor is attributed to the initial WP. 

To make clear the origin of the additional oscillatory contributions, we calculate, in 
sec. IIII[ the oscillation formula for a charged spin particle in the Sakata-Taketani Hamil- 
tonian formalism j^j], which is equivalent to the Klein-Gordon scalar wave equation. (The 



explicit analysis with mixed Klein-Gordon equation is made in Ref. |2Jj paying special at 
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tention to the relativistic initial value problem.) The oscillation formula in this case also 
possesses the additional interference terms between positive and negative frequency parts, 
very similar to the one obtained in the spin 1/2 case. From this example we will see that 
these interference terms are inevitable from a relativistic classical field theory perspective 
where covariance and causality is required. It is not specially associated to the spin degree 
of freedom. 

Another particular ingredient of neutrino oscillations can be included naturally within 
Dirac theory: the left-handedness of neutrinos created and detected through weak inter- 
actions. This fact, for a Dirac neutrino [2^, implies an additional probability loss due to 
conversion of left-handed neutrinos into right-handed neutrinos, which is possible because 
chirality is no longer a constant of motion for massive Dirac particles 2^|. Although previ- 
ous calculations |23| have shown an approximate contribution to this effect, we calculate in 
sec. Ill Al the complete effect. 

Concerning type (A) approaches, specifically the EWP description, we are interested to 
analyze further how is the propagation of intermediate virtual neutrinos. The framework 
where the investigations on first quantized approaches are made here is based on the cal- 
culation of the evolution kernels for free theories in presence of mixing. This enable us to 
deduce general oscillation probabilities in which there is explicit decoupling from the oscil- 
lating part (where all the oscillation information rests) and the initial wave packet. Another 
advantage of doing the calculations this way is that it resembles the propagator methods 
in covariant perturbation theory, which EWP approaches are based on. The free evolution 
kernel for fermions have a close relationship with the Feynman propagator used in QFT. 
What is common to both is that both particle and antiparticle parts contribute to the evo- 
lution or propagation. The necessity of the negative frequency part in the free evolution 
kernel is required from completeness and causality arguments but it also leads to the in- 
terference of positive and negative frequencies in flavor oscillation, treated in sees. |H] and 
TTT1 Then the question also arises in EWP approaches: are there contributions from both 
particles and antiparticles in the propagation of virtual neutrinos? In a simple microscopic 
scattering process, this question is meaningless since virtual particles are usually off-shell 
particles and must naturally have both contributions. However in EWP approaches the 
neutrinos propagate through macroscopic distances and, indeed, it can be shown 19 



that the virtual neutrinos are on-shell particles. There is no discussion, though, about the 
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possibility of neutrino and antineutrino contributions to the process; both can be on-shell. 
This investigation is carried on in sec. IIVI calculating explicitly the amplitude of produc- 
tion/propagation/detection process in an EWP approach. 

As a last task, we develop a simple, type (B), second quantized description of flavor os- 
cillation in sec. IIV Al using the free second quantized spin 1/2 fermionic theory in presence 
of mixing. This treatment has some similarities with the BV formalism but it does not 
require the introduction of flavor Fock spaces and Bogoliubov transformations. It means 
that the Fock space considered will be the one spanned by the mass eigenstates. Within 
this formalism it will be shown that the additional rapid oscillation contributions calculated 
through first quantized approaches do not survive the second quantization since only su- 
perpositions of particles (antiparticles) are used as initial neutrino (antineutrino) "flavor" 
states. Moreover, this property is not satisfied in the BV approach because the BV flavor 
states are mixtures of particle and antiparticle components; this is the ingredient responsible 
for a different oscillation probability [3j. 



II. FLAVOR OSCILLATION FOR DIRAC FERMIONS 

It is well known that the Dirac equation can give a significantly good description of a Dirac 
fermion if its inherent localization is much bigger than its Compton wave length; usually this 
is associated with weak external fields. For example, the spectrum for the hydrogen atom can 
be obtained with the relativistic corrections included (fine structure) [28|, p. 72]. One of the 
terms responsible for fine structure, the Darwin term, can be interpreted as coming from the 
interference between positive and negative frequency parts ( zitterbewequnq] of the hydrogen 
eigenfunction in Dirac theory compared to the nonrelativistic theory 24] . On the other hand 
a situation where the theory fails to give a satisfactory physical description is exemplified 
by the Klein paradox [3, p. 62]: the transmission coefficient for a electron moving towards 
a step barrier becomes negative for certain barrier heights, exactly when the localization of 
the electron wave function inside the barrier is comparable with its Compton wave length. 

Bearing in mind that first quantized approaches may fail under certain conditions we will 
treat in this section the flavor oscillation problem using the free Dirac theory in presence of 
two families mixing. The extension to treat three families of neutrinos is straightforward. 
A matricial notation will be used throughout the article for the first quantized approaches 
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to express the mixing. 

In matricial notation the mixing relation between flavor wave functions ^J(x) = 
(■0j e (x), (x)) and mass wave functions ^m( x ) = (V'i'C'Oj V'J ( x )) ^ s 

(cos# sin# \ 
*m(x) • (1) 
- sin# cos# / 

Each mass wave function is defined as a four-component spinorial function ip n {x, t), n — 1, 2 
that satisfy the free Dirac equation 

iJ^ n (x,0 = ff^nOM) , n = l,2, (2) 

where the free Hamiltonian is the usual 

= -ia.V + (3m n , n = l,2. (3) 

We will work in the flavor diagonal basis. This choice defines the flavor basis vectors 
simply as 

*j = (i,o), *; = (o,i), (4) 

while the flavor projectors are obviously 

P Va = i> a ul. (5) 

Actually, as an abuse of notation, the equivalence U ~ U <S> Id is implicit , as well as, 
Pu a ~ Pi/ Q <E> Id! the symbol 1^ refers to the identity matrix in spinorial space. 

The total Hamiltonian governing the dynamics of \l/ m is H D = diag(if-P, -f^f). From the 
considerations above, \P/(x, t) satisfy the equation 

i^fat) = UH D U- x *fat) . (6) 

The solution to the equation above can be written in terms of a flavor evolution operator 
K D as 

*/(x, t) = K D (t)^ f (x, 0) = J d 3 x' K D (x - x'; t)tf /(x # , 0) , (7) 



where 



i^(x - x'; t) = fl">(p; f) e^(—') . (8) 
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We can calculate K D (t) in any representation (momentum or position) as 

K D {t) = Ue- iuDt U- x 

cos 2 9 e- lH ? 1 + sin 2 6> e"*^ 4 — cos# smB(e~ lH ? 1 — e- lH ° l \ \ 

■ O D • ( 9 ) 

- cos6> sin6'(e" iH i 3 * - e -1 ^*) sin 2 # e ~ ljff f * + cos 2 ^ e" ,fl a f / 
The conversion probability is then 

Piy^v^t) = Jdx& f (x,0)K D \t)P Uli K D (t)y f (x,0) 

= Jdp^ e (p)(K^K^(p,t)^ e (p) , (10) 

satisfying the initial condition \I/J(x, 0) = (ip J e (x, 0) , 0) . Such initial condition implies, in 
terms of mass eigenfunctions, ^(x, 0) = ^(x, 0) = i[) Ve (x), as a requirement to obtain an 
initial wave function with definite flavor [9|. The function ?/v e (p) denotes the inverse Fourier 
transform of Vv e ( x ) (see Eqs. f)Al|) and (jA2)0 . 

Before obtaining the conversion probability for Dirac fermions, let us replace the spinorial 
functions */> n ( x ) by spinless one-component wave functions y?n(x) in the flavor wave function 
\I/J(x) — > (ip Ue (-x), (p u ^(x.)) and mass wave function ^m( x ) ~~ * (v^i( x ), ^(x)). We also replace 
the Dirac Hamiltonian in momentum space H E (p) (jSJ) by the relativistic energy E n (p) = 



p 2 + m 2 . Inserting these replacements into Eq. ()1U|) we can recover the usual oscillation 



probability [9 



23] 



&>{v e ^^;t) = J dx|z>Jtt /(x,t)| 2 

dp|^ e (p,*)vUp)P 



= lrfp^(p,t)l^(p)l 2 , (11) 

where ^/(x,0) T = (<^ e (x) T ,0), K% e {p,t) ee {K s ) 21 = - sin 9 cos 0( e ~ iE ^P) t - e - iE ^ 1 ) and 

&>(p, t) = sin 2 2fl sin 2 (A£( P )t/2) (12) 

is just the standard oscillation formula. The conversion probability (fTTj) in this case is then 
the standard oscillation probability smeared out by the initial momentum distribution. If 
the substitution |y5,, e (p)| 2 — > <5 3 (p — Po) is made the standard oscillation formula is recovered: 
it corresponds to the plane- wave limit. 

After we have checked the standard oscillation formula can be recovered for spinless 
particles restricted to positive energies in the plane- wave limit, we can return to the case 



of Dirac fermions. We can obtain explicitly the terms of the mixed evolution kernel Q by 
using the property of the Dirac Hamiltonian in momentum space H® 2 = E%(p) Id, which 
leads 

{KffKZfat) = &(p,t)[l-f(p)]l D 

+ sm 2 29f(p)sm 2 (Et)l D , (13) 

where 

/(P) = i[l - ^f^l , (M) 

and «^(p, t) is the standard conversion probability function f|12j) . A unique implication of 
Eq. (|13|). which is proportional to the identity matrix in spinorial space, is that the conversion 
probability (fTUj) does not depend on the spinorial structure of the initial flavor wave function 
but only on its momentum density as 

&>(v e ^;t) = Jdp{&(p,t)[l-f(p)] 

+ sin 2 2#/(p) sin 2 (^)}^t e (p)^ e ( p ) . (15) 

(The tilde will denote the inverse Fourier transformed function throughout this paper.) 
Furthermore, the modifications in Eq. (|15|) compared to the scalar conversion probability 
(jlljl are exactly the same modifications found in Ref. j3| and Ref. j3] after smearing out 
through a specific gaussian wave packet. 
The conservation of total probability 

^(i/ e ->i/ M ;t) + &>(v e ->v e ;t) = 1 , (16) 

is automatic in virtue of 

K£(p,t)K°(p,t) + K$(p 7 t)K°(p,t) = 1 D , (17) 

and the survival and conversion probability for an initial muon neutrino are identical to the 
probabilities for an initial electron neutrino because of the relations 

KflMK&fat) = K2(p,t)KZ(p,t) , (18) 
K£(p,t)K°(p,t) = K$(p,t)K°(p,t) . (19) 
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To explain the origin of the additional oscillatory terms in Eq. (|15p it is instructive to 
rewrite the free Dirac time evolution operator, in momentum space, in the form 

e -iHSt = e ~iE n t A D + + e iE n t A D_ ^ (2Q) 

where 

A° ± = \(ln±^-), (21) 

are the projector operators to positive (+) or negative (-) energy eigenstates of H®. By 
using the decomposition above (ffilj) . we can analyze in Eq. which contains the 
terms 

e iH?t e -iH?t = e iAEt A D +A D + + e -iAEt A D_ A D_ 

+ e i2 ^A? + Af_ + e" i2 ^Af_Af + , (22) 

plus its hermitian conjugate. Since A^A^- ^ 0, it can be seen that the rapid oscillating 
terms come from the interference between, e.g., the positive frequencies of the Hamiltonian 
Hi and negative energies of the Hamiltonian H®. One may think that by restricting the 
initial wave function to contain only positive energy contributions would eliminate the rapid 
oscillatory terms, as zitterbewegung disappears for superpositions of solely positive energy 
states in Dirac theory |28[, but it does not happen. The positive energy eigenfunctions 
with respect to a basis characterized by a mass mi necessarily have non-null components of 
negative energy with respect to another basis characterized by mi (this point is illustrated 
in appendix [Bj). Thus the rapid oscillatory contributions are an inevitable consequence 
of this framework and it is always present independently of the initial WP, if initially a 
definite flavor is chosen. However, its influence, quantified by the function f (p) in Eq. (fPfj) . 
is negligible for momentum distributions around ultra- relativistic values 22j. This rapid 
oscillatory terms will also be found for charged spin particle oscillations in the next section, 
with contributions slightly different from the ones obtained for spin 1/2 particles. 



A. Inclusion of Left-Handedness 



Until this point, we have been considering the oscillation of general flavor "particle num- 
ber" for general Dirac neutrinos. However, due to the left handed nature of weak interactions 
only left-handed components are produced and detected. To incorporate this fact into, for 
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example, the conversion probability in Eq. (fTUj) . it is sufficient to use initial left-handed WPs 
and replace the kernel of Eq. (fl"3"j) by the projected counterpart 

LK^(p,t)LK°(p,t)L = &> D (p,t)L 

1 .J„i '"I • ,- , "'2 



. sin 2 26 — ^ sin E x t - -==■ sin E 2 t L , (23) 
4 \Ei E 2 J 

where <^ D (p,t) = K$(p,t)K°(p,t) is the conversion kernel of Eq. ()13j) and L = (1 — 75) /2 
is the projector to left chirality. The conservation of total probability Ijlfij) no longer holds 
because there is a probability loss due to the undetected right handed component 

LK$RK£(p, t)L = \ sin 2 26 sin E x t - ^ sin E^tj \ , (24) 

where R — (1 + 75)/2 is the projector to right chirality. We can see that the probability loss 
(124)1 is proportional to the ratio m 2 /^ 2 which is negligible for ultra-relativistic neutrinos. 
The total probability loss for an initial left-handed electron neutrino turning into right- 
handed neutrinos, irrespective of the final flavor, is given by the kernel 

2 / m r. \ 2 



LK^RK^p,t)L + LK^RKg(p,t)L=[cos^6Q) sin 2 E x t + sin 2 ^) sin 2 E 2 t 



L . 

(25) 

To obtain the unphysical complementary kernels responsible for the conversion of right- 
handed component to right-handed and left-handed components, it is enough to make the 
substitution L <-> R in all formulas. 



III. FLAVOR OSCILLATION FOR SPIN 

The derivation of the usual conversion probability (jll)) takes into account only the positive 
frequency contributions. The mass wave function used to obtain Eq. (|TT|) corresponds to the 
solutions of the wave equation 

d 

1— p(x, t) = V-V 2 + mV(x, t) , (26) 

which is equivalent to the Dirac equation in the Foldy-Wouthuysen representation [29], 
restricted to positive energies. The evolution kernel for this equation is not satisfactory 
from the point of view of causality [30, p. 18], i.e, the kernel is not null for spacelike intervals. 



Moreover, the eigenfunctions restricted to one sign of energy do not form a complete set 
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24] 



To recover a causal propagation in the spin case, the Klein-Gordon wave equation must 
be considered. In the first quantized version, the spectrum of the solutions have positive 
and negative energy as in the Dirac case. However, to take advantage of the Hamiltonian 
formalism used so far, it is more convenient do work in the Sakata-Taketani (ST) Hamiltonian 



formalism 



2J] where each mass wave function is formed by two components 



<I>„(x./) - | ^ (X ' t} j . „: 1.2. 127) 

The components <p and \ are combinations of the usual scalar Klein-Gordon wave function 
4>(x) and its time derivative do4>(x). This is necessary since the Klein-Gordon equation is 
a second order differential equation in time and the knowledge of the function and its time 
derivative is necessary to completely define the time evolution. 

The time evolution in this formalism is governed by the Hamiltonian j^J 

H n T = -fa + iT 2 )^-+ml, (28) 

which satisfies the condition (H^ T ) 2 = (— V 2 + m^)l ST , like the Dirac Hamiltonian (|3|). The 
Tk represents the usual Pauli matrices and t ST is the identity matrix. 
A charge density jsjj can be defined as 

<5 n $ n = $tr 3 $ n = \<p n \ 2 - \ X n\ 2 , (29) 

which is equivalent to the one found in Klein-Gordon notation i(f>*do(j). Needless to say, this 
density (f2HJ) is only non-null for complex (charged) wave functions. The charge density 
is the equivalent of fermion probability density in the Dirac case, although the former 
is not positive definite as the latter. The adjoint $ = were defined to make explicit the 
(non positive definite) norm structure of the conserved charge 

J dxf> n (x, *)$ n (x, t) ee ($ n ,$„) = time independent . (30) 

Consequently, the adjoint of any operator Q can be defined as Cl = satisfying 
(fi$,$) = ($,fi$). Within this notation, the Hamiltonians of Eq. (j^HJ) is self-adjoint, 
H^ T = H^ T , and the time invariance of Eq. (|3*0*jl is assured. 

We can assemble, as in the previous section, the mass wave functions into = ($]", $J) 
and the flavor wave functions into fyj = $J ), satisfying the mixing relation *&f = U^ m . 
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The equivalence of U ~ U ® 1 ST and of P Ua ~ P Ua <S> 1 ST are implicit without modification 
in the notations. Then, the time evolution of can be given through a time evolution 
operator K ST acting in the same form as in Eq. (JZj). In complete analogy to the calculations 
from Eq. (JHJ) to Eq. (fTT)|) . we can define the conversion probability as 

= JdpIjp)K^K^(p,t)^ e (p) , (31) 

where \I//(x, 0) T = (<3> e (x) T ,0). The adjoint operation were also extended to ^ / = ^(1$ <E> 
r 3 ), where is the identity in mixing space. 

The information of time evolution, hence oscillation, is all encoded in 



KgKgfat) = ^(p,t)[l-/0*p)]lsT 

+ sin 2 2#/(/ip)sin 2 (£t)l ST , (32) 



where the function /(p) were already defined in Eq. (fTf]) and 



^JhUL + !?>). (33) 
V 2 m 2 mi 

The factor /i > 1 determines the difference with the Dirac case in Eq. (fT3*|) . The equality 
fi = 1 holds when mi = m-2, i.e., when there is no oscillation. 

IV. CONNECTION WITH QUANTUM FIELD THEORY 

The main improvement of the covariant approaches developed in sees. |H] and IIHI is that 
the propagation kernels governed by Dirac and Sakata-Taketani Hamiltonians are causal, 
i.e., are null for spacelike separations (see Eqs. ()A18)) and ()A19j) and Refs. Q,Q, Q)- On 
the contrary, the kernel of spinless particles restricted only to positive energies is not null for 
spaced 0.*o m the poiDtof view of relativistic classical field theories, a causal 

kernel guarantees, by the Cauchy theorem, the causal connection between the wave-function 
in two spacelike surfaces at different times j^]. 

To compare the IWP and EWP approaches it is useful to rewrite the Dirac evolution 
kernel for a fermion of mass m n , present in Eq. in the form [3, p. 89] 

K°{x-y) = £ [^[*;pK(!/;p) + <(x;pK(!/;p)]7o 

s J 
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= iS(x -y;m n )j , n = l,2, (34) 

where (x—y)° = t, (x—y)' 1 = (x— x') 1 when compared to the notation of Eq. (J2J). The spinorial 
functions u, v, are the free solutions of the Dirac equation and they are explicitly defined in 
appendix El (More familiar forms for the function S are also shown in appendix El) Clearly 
the function iS(x — y; m n ) = (0\{if) n (x) , ij) n (y)}\0) satisfies the homogeneous Dirac equation 
with mass m n (j2J and it is known to be null for spacelike intervals (x — y) 2 < 0, 

In contrast, the Feynman propagator iSp{x — y) appears in QFT. It is a Green function 
for the inhomogeneous Dirac equation obeying particular boundary conditions. The EWP 
approaches use this Green function for the propagation of virtual neutrinos. To directly 
compare the Feynman propagator to the kernel in Eq. (|3*3j) we can write iSf in the form 

iS F (x - y;m n ) = (0\T(^ n (x),ij n (y))\0) 

f Sp 

= E /^Kfop)<(2/;p)0(zo-2/o) 

s J 

-v s n (x;p)v s n (y;p)6(y -x )) . (35) 

Although the function Sp is called causal propagator, it is not null for spacelike intervals, and 
it naturally arises in QFT when interactions are present and treated in a covariant fashion. 
Equation (jHHj) shows that the propagator Sf describes positive energy states propagating 
forward in time and negative energy states propagating backward in time |2Ja, p. 91]. At 
a first glance, both neutrino and antineutrino parts of Eq. (|3*5|) seem to contribute to the 
space-time integrations present in covariant perturbation theory, as neutrino-antineutrino 
contributions in Eq. ()34j) have led to Eq. ([15)1 . 

In the following we will show in an EWP approach that for large separations between pro- 
duction and detection both neutrino and antineutrino parts may contribute as intermediate 
neutrinos for certain situations. 

We will follow the calculations made in Ref. j^J, using, instead of the scalar interaction, 
the effective charged-current weak lagrangian 

N=3 
i,a=l 

+v i (x)U* i ^Ll a (x)Jl(x)} (36) 
= ^ + &l , (37) 
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where G = \/2Gf and J M is the sum of any effective leptonic or hadronic current. The 
lagrangian is written only in terms of physical mass eigenstate fields, which coincides 



with flavor eigenstate fields on 
Suppose the process 



y for the charged leptons: l\(x) = e(x), h(x) = fi(x), .... 
33j where a charged lepton l a hit a nucleus A turning it into 
another nucleus A' with emission of a neutrino (this process happens around xa)- Sub- 
sequently the neutrino travels a long distance and hit a nucleus B which transforms into 
B' emitting a lepton lp (this process happens around xb)- The whole process looks like 
l a + A + B — > lp + A' + B' with transition amplitude given by 

(A'(p' A ),B'(p' B ),lp(pp)\S\A,BJ a ) . (38) 



The final states are momentum eigenstates while the initial states are localized [20( . The 
lowest order nonzero contribution of the scattering matrix S to Eq. (j3*B*jl is second order in 
the lagrangian (jSHJ). More explicitly, the term that contributes to the amplitude (J3~%|) comes 
from 

S (2) = t T ^ w f = -It{Z x + (39) 
~ -T(^)(^) (40) 
~ -G 2 I rf 4 xA^^(i,!/) , (41) 

J /3a 

where ( ) stands for space-time integration and 

JZ>p a (x,y) = ^^W^lf^iS^-^ijO^^!/)^!/): • (42) 

i 

In Eq. ()40p we kept only the mixed product and in Eq. ()41|) we kept from all possible terms 
in Wick expansion j3, p-180] only the term responsible for the transition of interest. 
Then the transition amplitude (J3%|) can be calculated as 

-G- 2 (A'(p' A ),B'(p' B ),lp(pp)\S^\A,B,l a ) 

d'y d'x (B'(p' B ) | J,{y)\B) (A'(p' A ) \ J\ (x)\A) 
x up{y,pp)^LY,UpiKi^F{y - x-m^L ^\l a {x)\l a ) (43) 

i 

J2 U PiKiAi- (44) 
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The initial states must be chosen in such a way that A, l a are localized around xa = if a-, x a) 
and B is localized around x B = {t B ,*- B ), since we are ultimately interested in large sepa- 
rations |x# — xa|. We can implement explicitly those localization conditions into the wave 
packets 

(B'(p' B )\J,(y)\B) = * / ie ip ^j; B '(q B)P ' B )fe(q B )e-^'fe^' (45) 

(2tc) 2 J 

(A'(p' A )\Jl(x)\A) = -Ly / ^e^-V / f A '(q A ,p^)^(qA)e-^- ( ™ ) (46) 
(0|Z„(x)|O = - ——g- / d^M^)e- iq ^ x - XA) , (47) 

(27r)2 J 

where ^ = dq{2E^))^l\ J BB '(q B ,p' B ) = (B'(p' B )\J^0)\B(q B )) and J^W Pa) = 
(A\p' A )\4(0)\A^A)). 

Following the calculations from Eq. (|43jl with the localization aspects of Eqs. (|45|) - (|4"7jl 
included, we arrive at 

A = J dq B J dq A J dq a J BB \q B ,p B )i) B (q B )tf A ' (cia,Pa)?Pa(cIa) 

e iq B .x Be i( qA +q a ).XA Uf} ( pf} ^ L J d 4 x d 4 y e iK e- y e- iKa - x iSi(y - x) 7 ^ Q (q Q ) > (48) 

where Kp = (ftp, Kp), n a = K a ) and 

Kp = pp + p' B - q B , K°p = Ep(pp) + E B ,(p' B ) - E b {c{b) , (49) 
K a = Pa ~ Pa + Qa , k° = E a (q Q ) - E A '(p'a) + £U(<U) • (50) 

By using the results of Eqs. (|C1J) and (|C2J) the expression between square brackets in 
Eq. fllBJ) gives 

J Anr 

x [u^k^u^k^d^i - mi) - Vi(-kur)vi(-kur)9(-LUi - m 4 )] , (51) 



where r = |y — x|, f = (y — x)/r, cUj = = k° and A;^ = y/uf — mf . The crucial point here 
is that, depending on the masses and momenta of the incoming particles, both neutrinos 
(mm) and antineutrinos (vv) can contribute to the amplitude (|4*H|) depending on the sign 
of its energy u>i, restricted to \u>i\ > m^ the off-shell contributions for Ui e [-mj, m,] are 
exponentially decreasing and then negligible for large distances (see appendix lUj) . We will 
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see in the following that antineutrino contributions in this case is possible and it corresponds 
to unphysical contributions. 

We are interested in large production-detection separations. It permits us to approxi- 
mate, as in Ref. [2Qj, r « _R + R.(y — x#) — R.(x — xa) and r ~ R, where R = |x# — x^| and 
R = (xb— x.a)/R. Such approximations inserted in Eq. (|51|) lead to momentum conservation 
on xa and xb vertices: 

AnR 

x [u^k^KjUi^KjO^uJi - rrii) - fc w R)vi(-fc w R)0(— a>i - m*)] . (52) 

At this point we have all the information to analyze whether the antineutrino part of the 
propagator contributes to the overall process. Neither of the isolated processes A + l a — > A' + 
i>i and B+i?i — > B'+lp are allowed if we calculate the transition amplitude for them separately 
using the weak Lagrangian (|36|). (For Majorana neutrinos they are strongly suppressed 
by helicity mismatch.) So far four- momentum conservation in both xa and xb vertices 
were automatically required from the calculations; among them the requirement of energy 
conservation for intermediate neutrinos with respect to the accompanying particles in vertex 
xa (&i = k°) and in vertex xb (^j = «$), is already implicit. The remaining are explicit 
in the delta functions of Eq. (|52j) . The on-shell condition for neutrinos (\uJi\ 2 — — mf) 
for long distance propagation was also automatic. What the calculations did not required 
is a definite sign for u>i, for all possible momenta constrained by the mentioned energy- 
momentum conservations. To analyze if and under what conditions both signs are possible 
is equivalent to study the kinematics of two-body to two-body scattering allowing the sign 
of one particle energy to be free. Putting in equations, for vertex xa, assuming the particle 
A at rest, we obtain from (pa — Pi) 2 = (pa' — Pa) 2 the neutrino energy 




The minimum value of right-hand side of Eq. (|53|) corresponds to last two terms equal to 
2m a MA'i which gives for the minimum 

min(^) = ^fWl ~ ( M a> ~ ™ a ) 2 + m 2 ] . (54) 
The values = K° a are bounded from below by the value in Eq. Imposing min(_Ej) > mi 
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and min(E'j) < — is respectively equivalent to 

M A - M A , > -(m a - mi) (55) 
M A -M A , < ~{m a + mi ) , (56) 

for M A > rrii and M A > > m a . It is clear that depending on the value of the masses, 
condition (|5Tj|) may be satisfied leading to antineutrino contributions to Eq. (jlH]) for a range of 
possible incoming momenta. Of course the condition is sufficient to exclude antineutrino 
contributions but it also excludes the cases where a threshold energy is required for the lepton 
I a to initiate the production reaction. Thus to prevent antineutrino contributions, it is better 
to adopt the weaker condition of restricting the sign of the energy of intermediate neutrinos 
Ui to be positive, keeping only the first term in Eq. (J51|) . Analogous analysis lead to possible 
momenta and mass values that allow k° < —rrii for vertex xb, still compatible with K° a = k%. 
Notice that condition K° a > rrii is exactly the kinematical condition to allow the production 
of physical neutrinos in x A and k ^ > rrii allow only the contribution of neutrinos with energy 
above threshold to trigger the detection reaction. The violation of these conditions implies 
in kinematically impossible contributions in production or detection. 

Restricted to condition > we can insert the expression above into Eq. (pIHj) which 
yields 

A = J S^2tt5(4 - K° a )6{ui - mi )^-e ik - R ~ i ^ t B-tA) e i{ Pk +P' B ).x B ^p' A .x A 
xu k {^p)^Lui{k^A)ui{k u A)YL^ a {^ a ) 



V T BB '(n W \ jMS^L jAA' (n n > \ M<Ia) 



\J Eb(s\.b) " ^E a (c\ a ) 



qs = Pg + p' b - k^R 

<1A = Pa _ lot + fc u R 



• (57) 



Notice that the step function 6{uJi — rrii) prevents non-physical neutrinos to contribute to 
the process. 

Particularly, if we use a unidimensional wave packet for the incoming lepton l a 

tpM = ^a(Qx, q y , Qz) = 5(q x )5(q y )ijj az (q z ) , (58) 
we obtain an amplitude analogous to Ref. j^J 



xUpiUpipp^LuiikuR) U^Uiik^Y^^ip 



APk+P' B )- x B p ip' A - x A 



<la=PaZ 
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v jBB' ( , \ ?MqB) jAA' (n W \ M<lA) 

s/E b [x\b) ^E a {<\a) 



. (59) 



where p a is the root of /(|q a | = p a ) = k° — k° = 0, which comes from energy conservation 
from the whole process; if there is no root the process is kinematically forbidden. The 
detection probability is proportional to the square of the amplitude (jHH)) integrated over the 
final phase space dp A dp B dp/3[2EA'(p , A)'^EB'(p B )2E/3(p0)]'' 1 . In particular, since P/3,PaiPb 
are fixed, the phases that differ for different intermediate neutrinos Ui are only k u R — C0i(t A — 
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ts) which is the same result obtained in Ref. 
velocity of particles A and B). 

So far we have shown in an EWP approach both processes in x A and xb should be 
considered real scattering processes with real neutrinos involved. The off-shell contributions 
are negligible to large distances and antineutrino contributions were explicitly excluded by 
eliminating the second term of Eq. (j51|) . These informations permit us to rewrite Eq. ()57|) 
in a slightly different form 

J d 4 x6(y - x)(A'(p' A ),is i {p)\JZ>l(x)e iP - XA \A,l a ) ,(60) 

where P = (H, P) is the energy- momentum operator. A change of notation were made here: 
in Eq. (jSHJ) the states \B) and \A, l a ) are centered around the origin while in Eqs. (jHBJ)-(j3Z!) 
they are respectively centered around xb and xa] the translation is explicitly performed by 
the translation operator e lRx . Additionally, the step function 8(y — x) is necessary to ensure 
that the contributions of points y around xb should always be after the contributions of 
points x around xa- By following the same steps from Eq. (|4H|) to Eq. (|58j) we can arrive 
from Eq. JSDJ to Eq. Ijgjty . 

Equation (J5UJ) shows us the amplitude of the entire process from production to detection 
in "decomposed" form (apart from the step function in time): the amplitude of production 
process multiplied by the amplitude of detection process summed over all possible inter- 
mediate real neutrinos of different masses m, and momentum p. (The sum over spins are 
implicit.) 
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A. A simple second quantized formulation 

Considering that only real neutrinos or antineutrinos (one of them exclusively) should 
travel from production to detection, the possibility to use the free second quantized the- 
ory for spin 1/2 fermions to describe flavor oscillations is investigated in this section. This 
simple, type B and QFT description of flavor oscillation phenomena guarantees only par- 
ticle or antiparticle propagation, keeping the nice property of giving normalized oscillation 
probabilities, like the first quantized examples treated in sees. |H] and IIHI 

To accomplish the task of calculating oscillation probabilities in QFT we have to define 
the neutrino states that are produced and detected through weak interactions. Firstly, we 
define the shorthand for the combination of fields appearing in the weak effective charged- 
current lagrangian (J5fijl 

v a (x) = U ai Vi(x) , a = e,/i. (61) 

We will restrict the problem to two flavor families and use the matrix U as the same in 
Eq. (JJJ. The mass eigenfields Vi(x), % = 1,2, are the physical fields for which the mass 
eigenstates |vj(p)) are well defined asymptotic states. The free fields Ui(x) can be expanded 
in terms of creation and annihilation operators (see appendix EJ and the projection to the 
one-particle space defines the mass wave function 

^(x; 5i ) = (0|^(x)|^: ft ) = E /dp4=<(x;p) , « = 1,2, (62) 

s J v ibji 

where 

Wi'9i) = ^J d P9 S i{p)Wi{P,s)) • (63) 

Since the creation operators for neutrinos (antineutrinos) can be written in terms of 
the free fields Pi(x) we can define the "flavor" states as the superpositions of mass 

eigenstates 

W-{g}) = u*M-9i) 

\vu'-{g}) = U ai \vi:gi) . (64) 

The details of creation are encoded in the functions 
We can also define 



*l>v a uefa {g}) = (Q\ve(x)\v a :{g}) 
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= U ei U* ai i) Vi {x\9i) , (65) 

where ip Ui ( x ) are then mass wave functions defined in Eq. (|62j) . We can see from Eq. 
that if ijj Ul (-x,t) = Vv 2 (x, t) = tp(x), for a given time £, ip UeVe (x,t) = ^(x) and Vv^ e (x, t) = 
due to the unitarity of the mixing matrix. 

Although this approach does not rely on flavor Fock spaces and Bogoliubov transfor- 
mations, we can use the same observables used by Blasone and Vitiello to quantify flavor 
oscillation j^: the flavor charges, which are defined as 

Qa{t) = j dx: z4(x,t)z/ Q (x,t) : , a = e,/i, (66) 

n 

where : : denotes normal ordering. Note that the Q e (t) + Q^{t) = Q is conserved [21], the two 
flavor charges are compatible for equal times, i.e., [Q e (t), Q^if)] = 0, and (u: {g}\Q\v- {g}) = 
±(v: {g}\v: {g}) for any particle state (+) or antiparticle state (-). Notice that in the second 
quantized version the charges can acquire negative values, despite the fermion probability 
density in first quantization is a positive definite quantity. The conservation of total charge 
guarantees the conservation of total probability (jl6j) . 

We can further split the flavor charges into left-handed (-) and right-handed (+) parts 

Q ( a\t) = y"dx:i/t(x,t)|(l±7 5 K(x,i) : , a = e,/i, (67) 

where + Qa = Q a - These components will be used to calculate the left-handed to 
right-handed transition. 

With the flavor charges defined, we can calculate, for example, the conversion probability 

^(v^v^t) = {v e :{g)\Q il {t)\v e :{g)) (68) 
= U^U^UiJdpe^-B^fag^togt) , (69) 

where the neutrino wave functions ip Vi are defined in terms of the function gt(p) in Eq. f)62|) . 
If we could equate the two mass wavefunctions in momentum space tp ui (p; gi) = Vv 2 (p; #2) = 
VVe(p) we would obtain, from Eq. ()69p. the standard two family conversion probability (fTTj) 

[Ve^Vrft) = J rfp^(p,t)^ e (p)^ e ( P ) , (70) 

where £P was defined in Eq. (I12j) . However, the equality can not hold as proved in ap- 
pendix (0): two wavefunctions with only positive energy components with respect to two 
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bases characterized by different masses can not be equal. Then, it is not possible to impose 
a flavor definite condition. Instead, we can write 

ft(p,«)=-^y*(p), (71) 

where ^(p) is the initial wave function associated to the neutrino of mass at creation, 
taking care to maintain the normalization J dp\gi(p)\ 2 = 1; any transition amplitude can 
be written in the form Eq. ([71)1 . In general ipi{p) = tp(p,rrii), and then, for small mass 
differences, 

Mp) ~ ^(P,m) ± —q- *l>(p,™>) , (72) 

where m = (mi + rri2)/2 and Am = ni2 — m\. Keeping only the first term, ip(p, m) = ^(p), 
we obtain from Eq. (jEHJ), 

&{u e -+vrit) = jdp^{pM\p)[t-\K x ^p)-\K^{p)]^{p) 

/~ /\ 77? ~ 

rf P V t (p)[/(p)cos(AEt) -i__ 7 .p S in(AEt)]^(p) (73) 

Notice that in this case, the conversion probability is non-null for t = 0, 

^ e -^;0) = isin 2 2fl fdpf{p)ft{p)iP{p) , (74) 

which imply a direct lepton flavor violation in creation. However, since /(p) ~ (Am) 2 / (4p 2 ) 
for ultra- relativistic momenta, the violation is hopelessly feeble for direct measurement. 
Among the deviations of the conversion probability (J7H|) compared to the standard one (JZOJ), 
only the last term is of order Am/E, the rest is of order (Am/E) 2 (the contributions of Aj_ 
can be estimated by [vl(— p, s)«i(p, s')} 2 ~ p 2 [Am + A£'] 2 /[(m 1 + E\)(mx + Ex)]). Even so, 
Am/E ~ 10~ 9 for Am 2 ~ 10~ 3 eV 2 , m ~ l/2eV and E ~ lMeV, which can not be seen in 
actual oscillation experiments. Nevertheless, it is important to note that the knowledge of 
Am in conjunction with Am 2 gives information about the absolute mass scale because of 
Am 2 = 2mAm. 

Using Qo - -* of Eq. (|57|) instead of Q a in Eq. and t/>(p) = L^(p) in Eq. (|70|) we obtain 



^ez-^ ; o = /"dp \^S-np, t) + \ sm 2 M 11 ^-^ 

v M ' 7 7 H \ kE x E 2 v 74 V2E! 2£ 2 



^(pMp) • (75) 
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The total probability loss from the conversion of initial left-handed electron neutrino to 
right-handed neutrinos yields 



dp 



cos 2 6 



mi 



+ sin 2 6» 



m 2 



^(pMp) • (76) 



- 2Ei / \ 2E2 - 

Notice Eq. (fTtjj) does not depend on time in contrast to its first quantized analog in Eq. (|23j) . 
Other conversion and survival probabilities can be obtained from Eq. (JT7jl and 



(77) 



The exchange of L <-> R does not modify the formulas, provided we also change the chirality 
of the initial wave function. 

For completeness we calculate the additional conversion probabilities including the second 
term of Eq. (I72|) 



6&>( 



, 2tf— /dp^(p) 



Ho 



Hi 



+ 



Am 



+(Ai+ + A 2+ - f(p))i sin AEt 



IE2 2Ei IE1E2 

^t^(p) + h - c 
am 



TP 



i sur > 2 0^ I dp ft (p) 



+ J z sin AEt 



— ^(P 

dm 



IE1E2 

which have terms of order Am and (Am) 2 . 

To calculate the conversion probability for antineutrinos u, 



2i?2 2£?i ' 
+ h.c. , 



(78) 



(79) 



Vp, it is enough to use 



(80) 



instead of Eq. (J7T|) . replace t — ► — t and V>(p) - > t^F ^(p) i n t ne expressions corresponding to 
neutrinos (jRHjl — (j75)l . or apply charge conjugation ftp) — > — 7oC^*(p). These prescriptions 
can be inferred from direct comparison to the calculations and beware that the definition of 
antineutrino states are defined with gf*(p) (|A29[) . 

The formulas obtained in this second quantized version does not have the interference 
terms between positive and negative energies like in Eq. (|15jl. Such interference terms are 
absent because the possible mixed terms like &2(p)o4(p)|0) are null for an initial "flavor" 
state superposition that contains only particle states (or only antiparticles states). Notice 
that the irrelevance of the initial spinorial structure no longer holds in this second quantized 
version, which can be seen, for example, in Eq. ([73)1 . 
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V. DISCUSSION AND CONCLUSIONS 



Using the Dirac equation which is a relativistic covariant equation we obtained oscillation 
probabilities respecting causal propagation. The oscillation formulas obtained had additional 
rapid oscillating terms depending on the frequency E\ + E%, with respect to the usual 
oscillation formulas with wave packets. Such additional oscillatorial character seemed to 
have its origin in the intrinsic spinorial character used. However, we have seen that such 
terms also appear in the charged spin particle oscillations. In fact the rapid oscillatorial 
terms arise from the interference of positive and negative frequency parts of the initial 
WP and they are always present independently of the initial wave packet if initially a 
flavor definite condition were imposed. In addition, within Dirac theory, we have shown the 
detailed spinorial character of the initial wave function was irrelevant for flavor oscillation, 
independently of the momenta involved. The inclusion of the left-handed nature of the 
created and detected neutrinos could also be simply achieved. It is important to stress 
that the modifications found in this context would have tiny observable effects to the flavor 
oscillation of ultrarelativistic neutrinos. 

Regarding second quantized approaches (sec. IIV|). in particular, EWP approaches, we can 
compare the propagators used in the latter to the evolution kernels in IWP approaches. Both 
the free evolution kernel and the Feynman propagator for fermions contain the contribution 
from particle and antiparticles. From this perspective, EWP approaches could also contain 
both contributions from neutrinos and antineutrinos, as in the first quantized approaches 
presented in sees. HT1 and IIHI To analyze this point, an EWP calculation were carried out 
explicitly in sec. HVI following Ref.l^. Restricted to a case where only neutrinos would be 
present, the calculation showed that the antineutrino contribution were not excluded auto- 
matically in the formalism but a subsidiary condition could be necessary: the sign of the 
frequency of the intermediate neutrinos should be restricted to be positive. In such case, 
there can be interference terms between positive and negative frequencies, possibly yielding 
rapid oscillation terms similar to the ones obtained in Dirac theory of sec. El However, it 
should be stressed that the origin of the interference between positive and negative terms are 
different in the first quantized Dirac theory treated in sec. |H] and in the EWP (second quan- 
tized) treated in this sec. IIV1 The former comes from causality and completeness arguments 
in a classical field theory perspective, while the latter have its origin in the consideration 
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of non-physical contributions. The restriction implied by the subsidiary condition automat- 
ically guarantee that: (i) only real neutrinos that are kinematically allowed in production 
contributes and (ii) in detection, only neutrinos with energies above threshold to trigger the 
detection reaction contribute. Otherwise, kinematically forbidden reactions in production 
or detection could be possible through exchange of virtual antineutrinos carrying negative 
energy. The overall energy-momentum conservation, though, is always respected (smeared 
out through finite momentum distributions) through production/propagation/detection pro- 
cesses. Since the presence of both neutrino-antineutrino contributions is common to all EWP 
approaches, the subsidiary condition necessary in the EWP approach analyzed is possibly 
necessary in any approach with virtual neutrino propagation. (Unless a stronger condi- 
tion like Eq. (|HH)l is already implicit.) For example, in Eq. (14) of Ref. Q, the subsidiary 
condition (for antineutrinos) is satisfied because the detection reaction is an elastic scatter- 
ing. (Although the detection electrons are off-shell (bound state), the subsidiary condition 
is valid.) An important remark in this respect is that the calculations of the production 
and detection amplitudes as separate processes take automatically into account the physi- 
cal kinematical conditions (i) and (ii) through the energy-momentum delta functions. It is 
also important to stress that the result obtained here did not depend on particular wave 
packets neither on the particular interaction used. The interesting point is that by impos- 
ing such subsidiary condition beforehand enables us to write the transition amplitude for 
the entire production/propagation/detection process in a "decomposed" form, with simple 
physical interpretation. A more detailed discussion about the decomposition of the process 
in separate production, propagation and detection processes can be found in Ref. |2jJ. Similar 
conclusions can be drawn for the case where only antineutrinos should propagate: the sign 
to be chosen should be negative. A realistic EWP approach for antineutrino propagation is 
given in Ref. [l9|. To conclude this part, EWP approaches for neutrino oscillations require 
for consistency, but do not automatically imply, real intermediate neutrinos or antineutrinos 
exclusively. 

All the properties discussed above about EWP approaches suggest that the description 
of two macroscopically distant scattering processes (production and detection) as a single 
scattering process described by a single scattering matrix have to be treated with care. 
We can be confident about the use of the S matrix to describe any microscopic event 
through perturbation theory to any order of expansion (if the theory is renormalizable) , 
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but the extension to macroscopically distant reactions is not automatic. Actually, if the two 
processes are indeed not causally connected it can be proved that the S matrix decomposes 
as the product of the two S matrices for the two distant and independent processes [35| . 

From the considerations above, the positive and negative interference terms in the first 
quantized approaches considered seem unphysical. To support this idea, it was shown in 
sec. HV Al that a simple second quantized, type B, and IWP-like, approach could be devised 
using the second quantized free theory maintaining the simple properties of IWP approaches 
but eliminating the undesirable interference of positive and negative frequencies that was in- 
evitable in the relativistic quantum mechanical context. On the other hand, new ingredients 
such as the direct flavor violation in creation and deviations from the standard oscillation 
formula were found. The deviations include the probability loss due to the conversion of 
left-handed neutrinos to right-handed neutrinos. Unfortunately, those new effects are tiny 
either because of the small mass difference or the ultrarelativistic nature of neutrinos and 
they are not feasible for direct observation in actual oscillation experiments. 
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APPENDIX A: NOTATION AND DEFINITIONS 

The (scalar, spinorial or ST) wave functions related by Fourier transforms are denoted as 

= /^(P)e ip ' x , (Al) 

#P) = JdxvMe-**. (A2) 

The tilde denotes the inverse Fourier transformed function. 

Using the property of the Dirac or ST Hamiltonian, H% = (p 2 + m n ) 2 l, we can write the 
evolution operator in the form 

e ~iH n t = cos (£ nt ) _ i_j± s in(£U) , (A3) 

E n 

where the momentum dependence have to be replaced by — zV in coordinate space. 
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The free neutrino field expansion used is (i = 1, 2) 

= eJt^- [<(z;pH(p, s) + <(x;p)6t(p, s)] , (A4) 

where the creation and annihilation operators satisfy the canonical anticommutation rela- 
tions 

{0i(p, r), a}(p', s)} = <%5 rs 2^(p)S 3 (p - p') , (A5) 
{fcfoO^p'.s)} = ^^(p^p-p') ; (A6) 

all other anticommutation relations are null. The functions u, v are defined as 

e — ipi-x 

ut(x;p) = <(p)^i7i , (A7) 
<(P) = /7? ■ — u , (A8) 

,/ \ m i -£ i 7° + P-7 

u (p) = 7= w , ( Al 0) 

+ 

where p^.x = Ei(p)t — p.x and they satisfy the properties 

u r u s = u^Uq = -v r v s = v$v% = 6 rs , (All) 
v$u s = v$v' = Vr.s, (A12) 
E«?(P)**(P) =^ + «» = 2^(p)Af + (p) 7 , (A13) 

s 

E«?(P)«<'(P) = ^-™« = 2^(p)Af_(-p) 7 • (A14) 
The Feynman propagator for fermions is 

iS F (x-y) = (0\T{1>{x)${y))\0) (A15) 
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- e"^"^ (A16) 
(2vr) 4 i> - m + ze 

= (z^ + m)iAi?(x — y; m) . (A17) 

The function S in Eq. (|34|) and its equivalent for the Sakata-Taketani Hamiltonian can 
be written as 

iS(x;m) = (ift + m)iA(x; m) , (A18) 
V 2 

K ST (x; m) = [id t - 7^( T 3 + ir 2) + m 2 }iA(x; m) , (A19) 



iA(x;m) = — — / <fp5{p 2 - m 2 )e(p )e ^ 
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f '' P [e- ip - x - e +ip - x ] . (A20) 



(2TifJ 2E p 
The free neutrino eigenstates are defined as 

h(p, S )> ee ^^|0) (A21) 

= J dxv}(x)\O)^0 , (A22) 

|^( P)S )) ee ^f=^0> (A23) 

= J dK^^-^xM , (A24) 

whose normalization is (^(p, r)|z/j(p', s)) = 5ij5 rs 5 3 (p — p'). The same normalization is valid 
for the antiparticle states. The states with finite momentum distributions are defined as 

W-9) = E Jdpg s (p)\v z (p,s)) (A25) 

= J dxv\(x)\W Vi {x) , (A26) 

Mx) ee £ JdP9 s (p)^^, (A27) 

e-^>,:y) = y rfxz/l(x,0)|0)^(x,t) , (A28) 

W-9) = E Jdpg s *(p)\Pi(p,s)) (A29) 

= y dx^^^^lO) , (A30) 

iM*) = E Jdpg s (p) V -^0t, (A31) 



^:<?) = y'dx4 i (x,t) I / i (x,0)|0) . (A32) 



APPENDIX B: DECOMPOSITION WITH RESPECT TO TWO BASES 

It is possible to decompose a given spinorial function ^( x ) i n terms of bases depending 
on different masses mi and m 2 . Equating 

^ (x) = Im [<(x; p) ^ (+)(p ' s) + v * (x; p) ^ ( " )(p ' s)] ' 4 = h 2 ' (B1) 

the expansion coefficients can be obtained 

^ (+) (p,s) = y'rfx< t (x;p)^(x) , (B2) 

gt\p,s) = y'rfx< t (x;p)^(x) . (B3) 
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From Eq. (jB3|) we see that imposing the conditions 



,4, 



pM-p) 
pM-p) 



for all p, lead to the equivalent conditions 



[(m l + E 2 )-(m 2 + E 2 )}v s ^(-p) 
1 1 



mi + E 2 m 2 + E 2 



^o f 7-P^(-p) 



0, 
0, 



, s 
, s 



1,2 
1,2 



(B4) 
(B5) 



(B6) 
(B7) 



where the property of Eq. (jA10|) and 70 Vq = —Vq were used. In case mi 7^ m 2 , we can use 



the decomposition ip(p) = ^+(p) +V L -(p)> where ^>±(p) 
Eqs. f|B6j) and (JB7j) the conditions 



^1 ± 7o)V ; (p)/2, and obtain from 



«oV-(-p) 



Uo t o--pV'+(-p) 








1,2, 
1,2, 



(B8) 
(B9) 



where the properties 7 = 7075c and u s Q = 75^ were used in Eq. (IB 9)1 . The equations (JBf 
and (|B9|) are only satisfied if VM-(p) = V'-(p) — 0, since cr.p has only non-null eigenvalues 
and it commutes with 70 . It is easier to reach this conclusion in the helicity basis {u , } 



characterized by cr.puQ 



±|p|tiQ , but the result is basis independent. 



APPENDIX C: INTEGRALS 

Splitting the Feynman propagator into positive and negative frequency parts iSp(x) = 
iS (+) (x) + iS ( ~\x) the following integrals give us 

dte iujt iS <+) (r } t;m) = - k^r.-y) + m}6(uj - m) (CI) 

dt^HS^Hr^m) = H)^[(w7o-/Ur. 7 ) +m]0(-u-m) , (C2) 



where k w = \J uj 1 — m 2 , r = rr, the conditions mr, k u r ^> 1 were assumed and terms behaving 
as 1/r 2 were neglected. 

To illustrate the calculations, Eq. (jCl|) is obtained by 



V ; (2tt) 3 J 2E(p) u-E(p)+ie V ; 

1 —i f°° p ( . E(p) r y + m 



(2vr)2 2r7_ 0O E(p) I ' E{p) - u - ie 

r sinpr tp(f.7) ( 
+ cospr h=rr^ f • C4 
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In Eq. ()C3|) the following identity is used 



o 



To get to the result of Eq. (jCl|) it is necessary to split the functions sinpr and cospr in 
Eq. (|C4j) into exponentials and, for the e tpr part, integrate along a closed path formed by 
a half semicircle in the upper-half complex plane going round the branching line [im, ioo) 
(for the e~ ipr part take the path reflected by the line defined by Rez = 0). The contribution 
from the paths beside the branching line yields a function which decreases more rapidly than 
e~ mr and it is negligible for mr ^> 1. The contributions for — m < u < m give a function 
with dependence e~^ r which is also negligible for large separations r. 
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